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ABSTRACT 

We compare quantitative estimates of the filamentarity of the galaxy distribution 
in seven nearly two dimensional sections from the survey against the predictions of 
ACDM N-body simulations. The filamentarity of the actual galaxy distribution is 
known to be luminosity dependent. It is also known that the filamentarity of the 
simulated galaxy distribution is highly sensitive to the bias, and the simulations are 
consistent with the data for only a narrow range of bias. We apply this feature to sev- 
eral volume limited subsamples with different luminosities to determine a luminosity 
bias relation. The relative bias b/b* as a function of the luminosity ratio L/L* is found 
to be well described by a straight line b/b* = A + B {L/L*) with A = 0.833 ± 0.009 
and B = 0.171 ± 0.009. Comparing with the earlier works all of whic h use ratios of 
the tw o-poi nt statistics, w e find that our results are consistent with iNorberg et ail 
(l200lh and iTegmark et all (|2004f ). while a steeper luminosity dependence found by 
iBenoist et al.l (|1996( ) is inconsistent. 



Key words: methods: numerical - galaxies: statistics - cosmology: theory - cosmology: 
large scale structure of universe 



1 INTRODUCTION 

Filaments are possibly the most prominent visible feature 
in galaxy redsh ift surveys. The Las Campanas Redshift 
Survey (LCRS; IShectman et~aT1 [l996) has been analyzed 
showing the observed filamentarity to be statistically sig- 
nificant on length-scales upto ~ 80/i ~ Mpc and not beyond 
(Bharadw ai. Bhavsar fc S hcth 2004]), establishing the fila- 
ments to be the largest known statistically significant coher- 
ent structur es in the galaxy d istribution. The Sloan Digital 
Sky Survey l|York et al.l[20 00) is currently the largest galaxy 
redshift survey. Studies using two nearly two-dimensional 
(2D) sections from the Sloan Digital Sky Sur vey Data Re- 
lease one (SP SS DR1: lAbazaiian et all hoO&D confirm the 
earlier result l|Pandev fc Bharadwail I l2005l ~ This analysis 
also shows the filamentarity to be luminosity dependent with 
the brighter galaxies having a more compact and less fila- 
mentary distri bution as compared to the fa int ones. A more 
recent study (IPandev fc Bharadwail 120061) using a larger 
data-set from the Sloan Digital Sky Survey Data Release 
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Four (SDSS DR4: lAdelman- McCarthy et al.ll2006l ) confirms 
the luminosity dependence and shows this effect to be con- 
siderably enhanced for galaxies brighter than L* . This study 
also shows the filamentarity to be dependent on the colour 
and morphology of the galaxies. 

The possibility that the galaxies are a biased tracer 
of the underlying dark matter distribution is now quite 
widely accepted. A study comparing the filamentarity ob- 
served in the LCRS with the prediction of ACDM dark mat- 
ter N-body simulations finds the N- body predictions to be 
very sensitive to the bias parameter (jBharadwai fc Pandevl 
|2004| ). Increasing the bias causes the simulated galaxies to 
get concentrated into smaller regions which correspond to 
the peaks in the dark matter distribution. The large-scale 
filamentarity of the simulated galaxy distribution is found 
to decrease with increasing bias. The filamentarity of the 
simulated galaxy distribution and the LCRS galaxies in the 
absolute magnitudes range —21.5 < M < —20.5 are found to 
be consistent for a bias b = 1.15. An unbiased ACDM model 
(6 = 1) and a large bias parameter (6 = 1.5) are both ruled 
out. Comparing the filamentarity in biased N-body simula- 
tions with the actual galaxy distribution provides a novel 
technique to determine the galaxy bias. 
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The large coverage and high galaxy number density of 
the SDSS permits us to study the filamentarity as a function 
of galaxy luminosity. In this Letter we compare the filamen- 
tarity of galaxy samples in different luminosity bins drawn 
from the SDSS DR4 against the predictions of biased N- 
body simulations and use this to determine the bias as a 
function of luminosity te. a luminosity bias relation. 

The luminosity dependence of galaxy clustering 
is an important issue with considerable implications 
for theories of galaxy formation. Observations show 
that the luminous galaxies exhibit a stron ger clus- 



tering than their fainter counterparts (e.g. Hamiltonl 
19981. iDavis et all Il988l, IWhite. Tullv fc Davisl 



Park, Vogclcv, Geller & Huc hral 1994. lLovedav et al 



Guzzo et al.lll997UBenoist et alj|l99rl iNorberg et allhoOll. 



1988 



1995, 



Zehavi et alj 2005h . It has remained difficult to quantify 



this effect through a luminosity bias relation because of the 
limited dynamical range of even the largest redshift surveys. 
Traditionally the luminosity bias relation has been deter- 
mined by comparing observations of the two-point statistics, 
namely the correlation function or the power spectrum in 
different luminosity bins. This allows the relative bias 6/6* 
to be studied as a function of the luminosity ratio L/L* 
where 6* is the bias corresponding to the galaxies with the 
characteristic luminosity L* or equivalently the characteris- 
tic mag nitude M* . To our k nowledge there are three earlier 
results. iBenoist et al.l (|l996l ) have an alyzed the luminosity 
dependence of the bias in the SSRS2 i|da Costa et al.111994 ^ 
which has been found to be well fitted b y a luminosity 
bias relation (6/6*) = 0.7 + 0.3 (L/L*) JPeacock et all 
20011). The an a lysis o f the 2dFGRS (jColless et al.ll200ll ) by 



Norberg et al.l (|20011 ) yields a more modest luminos i ty de - 
pendence (6/6*) = 0.85 + 0.15(L/L* ).lTegmark et all (|2004T ) 
have analyzed the SDSS (lAbazaiian et all 2004h to obtain a 
relation (6/6*) = 0.895+0.15(L/L*)-0.04(M-M*). Unlike 
these estimates, the present analysis uses a global property 
namely the filamentarity, to determine the luminosity bias 
relation. 

There currently exist several ways to directly determine 
the bias. Measuring the redsh ift space distortion p arame- 
ter (3 = Qt e /b (e.g 2dFGRS. llTawkins et al.ll2003l ; SDSS, 
iTeemark et al 1 l2004h in combination with an independent 
determination of Q mo allows 6 to be determined. The bis- 
pectrum (2dFGRS, IVerde et al.ll2002l ) provides a technique 
to determine the bias from redshift surveys without the 
need of inputs from other observations. A combination of 
we ak lensing a n d the SDSS galaxy survey has been used 
bv lSeliak et all (|2005f ) to determine the bias. We note that 
most of these methods require large galaxy samples with 
a high number density. It has not been possible to apply 
these techniques to subsamples with different luminosities 
and determine a luminosity bias relation. 



2 DATA AND METHOD OF ANALYSIS 

For the present analysis we use seven non-overlapping strips 
each with sky coverage 2° x 90° which lie entirely within 
the survey area of SDSS DR4. The se strips are identical 
in sky coverage as the ones used in iPandev fc Bharadwail 
(2006) and are shown in Figure 1 of that paper. For each 
strip we have extracted five different volume limited sub- 



samples, 'Bin 1' ... 'Bin 5', applying the absolute magnitude 
and redshift limits listed in Table 1. The thickness of the 
resulting subsamples increases with redshift. For our analy- 
sis we have considered a smaller region of uniform thickness 
corresponding to the value at the lowest redshift. For each 
bin the area and thickness are listed in Table 1. For all the 
subsamples, the thickness is much smaller than the other 
two dimensions and hence it is collapsed along the thickness 
resulting in a 2D distribution (Figure [TJ. For each luminos- 
ity bin the galaxy number density varies slightly across the 
seven strips and the average density along with the 1 — a 
variation is shown in Table 1. It may be noted that the bins 
span AM = 1 in absolute magnitude and are in order of in- 
creasing luminosity. The characteristic magnitude M* has a 
value M* = -20.44±0.01 for the SDSS (|Blanton et al.ll2003T l 
and the A/* galaxies appear in Bin 3 of our sub-samples. The 
galaxy number density falls considerably in the bins contain- 
ing galaxies brighter than M* . 

We simulate the dark matter distribution using a 
Particle-Mesh (PM) N-body code. The simulations use 256 3 
particles on a 512 3 mesh, and they have a comoving vol- 
ume [921.6/i _1 Mpc] 3 . We use (Cl m0 ,n A0 ,h) = (0.3,0.7,0.7) 
for the cosmological parameters along with a ACDM power 
spectrum wi th spectral index n s — 1 and normalization 
0.84 (ISpergel et al l 120031). 



<T8 

ing 



scheme 



A "sharp cutoff" bias- 



Cole. Hatton fc Weinbergl 1 199^ was used to 



extract particles from the N-body simulations. Identifying 
these particles as galaxies, we have galaxy distributions that 
are biased relative to the dark matter. The bias parameter 
6 of each simulated galaxy sample was estimated using the 
ratio 



6 = 



(i) 



where £ 9 (r) and £(r) are the galaxy and dark matter two- 
point correlation functions respectively. This ratio is found 
to be constant at length-scales r > 5/i _1 Mpc and we use the 
average value over 5 - 40/i _1 Mpc. We use this method to 
generate galaxy samples with bias values in the range 1 to 
2 in steps of 0.1. Galaxy distributions with bias < 6 < 1 
were generated by adding randomly distributed particles to 
the dark matter distribution. In addition to the bias, the 
filamentarity of the simulated galaxy distribution is also ex- 
pected t o depend on the de gree of non-linearity in the evolu- 
tion (see I Verde et al.ll2002l for effects on the bispectrum). To 
asess this effect we have run simulations using a% = 1 and 
the best fit cosmological p arameters determin ed from the 
SDSS+WMAP 3 year data (|Spergel et al.ll2006l . a 8 = 0.77), 
and find that variations in this range do not cause a statis- 
tically significant change in the filamentarity. 

The peculiar velocity effects were included to produce 
galaxy distributions in redshift space. We have used three 
independent realisations of the N-body simulations, and for 
each luminosity bin we have extracted six different simu- 
lated data-sets from each simulation. This gives us a total 
of eighteen simulated galaxy distribution for each luminos- 
ity bin. The simulated galaxy distributions (Figure [TJ have 
the same area, thickness and the mean number density as 
those of the corresponding bins (Table 1). The simulated 
data were analyzed in exactly the same way as the actual 
data. 

All the sub-samples that we have analyzed are nearly 
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Table 1. This shows the absolute magnitude and rcdshift limits for the different volume limited subsamplcs analyzed. The median 
magnitude, area, thickness and average galaxy number density with 1 — a variations from the 7 strips are also shown. 



Bin Magnitude range Rcdshift range Area Thickness Density 

10 4 /i- 2 Mpc 2 ft^Mpc 10" 2 /i 2 Mpc- 2 



Bin 1 


-18.5 > M r 


> 


-19.5 


0.022 < z < 0.061 


2.21 


2.33 


1.93 ± 0.16 


Bin 2 


-19 > M r 


> 


-20 


0.028 < z < 0.075 


3.28 


2.91 


2.33 ±0.2 


Bin 3 


-20 > M r 


> 


-21 


0.043 < 2 < 0.114 


7.41 


4.52 


1.77 ±0.15 


Bin 4 


-21 > Mr 


> 


-22 


0.067 < 2 < 0.169 


15.5 


6.94 


0.67 ± 0.04 


Bin 5 


-21.5 > M r 


> 


-22.5 


0.083 < 2 < 0.2 


20.9 


8.55 


0.25 ±0.01 
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Figure 1. This shows the galaxy distribution in one of the SDSS 
strip in Bin2 along with some of the simulated slices with bias 
values as indicated in the panels. 

two dimensional. The samples were all collapsed along the 
thickness (the smallest dimension) to produce 2D galaxy dis- 
tributions (FigureCD). W e use the 2D "Shapefinder" statistic 
ijBharadwai et al.|[200(il ) to quantify the average filamentar- 
ity of the patterns in the resulting galaxy distribution. A 
detailed discussion of t he me thod of analysis is presented in 
iPandev &; Bharadwail l|2006h (and references therein), and 
we highlig ht only the salient features here. The rea der is 
referred to ISahni. Sathvaprakash.fe Shandarinl fl998) for a 
discussion of Shapefinders in three dimensions. 

The galaxy distribution is represented as a set of Is on a 
2-D rectangular grid of spacing 1 h~ Mpc x 1 h~ 1 Mpc. The 
empty cells are assigned a value 0. We identify connected 
cells having a value 1 as clusters using the 'Friends-of- Friend' 
(FOF) algorithm. The filamentarity of each cluster is quan- 
tified using the Shapefinder T defined as 



(P 2 - 165) 



(2) 



(P-4Z) 2 

where P and S are respectively the perimeter and the area 
of the cluster, and I is the grid spacing. The Shapefinder T 
has values and 1 for a square and filament respectively, 
and it assumes intermediate values as a square is deformed 
to a filament. We use the average filamentarity 



F 2 = 



(3) 



to asses the overall filamentarity of the clusters in the galaxy 
distribution. 

A problem arises because the distribution of Is corre- 
sponding to the galaxies is rather sparse. Only ~ 1% of 



the cells contain galaxies and there are very few filled cells 
which are interconnected. As a consequence FOF fails to 
identify the large coherent structures which are visually dis- 
cernible as filaments in the galaxy distribution (Figure [Q. 
We overcome this problem by successively coarse-graining 
the galaxy distribution. In each iteration of coarse-graining 
all the empty cells adjacent to a filled cell are assigned a 
value 1. This causes clusters to grow, first because of the 
growth of individual filled cells, and then by the merger of 
adjacent clusters as they overlap. Coherent structures ex- 
tending across progressively larger length-scales are identi- 
fied in consecutive iterations of coarse-graining. So as not to 
restrict our analysis to an arbitrarily chosen level of coarse- 
graining, we study the average filamentarity after each iter- 
ation of coarse-graining. The filling factor FF quantifies the 
fraction of cells that are filled and its value increases from 
~ 0.01 and approaches 1 as the coarse-graining proceeds. We 
study the average filamentarity F2 as a function of the filling 
factor FF (Figure [3]) as a quantitative measure of the fila- 
mentarity at different levels of coarse-graining. The values 
of FF corresponding to a particular level of coarse-graining 
shows a slight variation across the data-sets. In order to 
combine and compare the results from different data-sets, 
for each data-set we have interpolated F2 to 7 values of FF 
at an uniform spacing of 0.1 over the interval 0.05 to 0.65. 
Further coarse-graining beyond FF ~ 0.65 washes away the 
filaments and hence we do not include this range for our 
analysis. 

For each luminosity bin (Table 1) the seven different 
strips were used to determine the mean and the variance- 
covariance matrix for F2 as a function of FF. This was com- 
pared against the mean and the variance-covariance matrix 
estimated from eighteen simulated data-sets for each value 
of the bias. Bias values in the range 0.5 — 2.0 were consid- 
ered. For each v alue of the bias w e use the Hotelling's two 
sample T 2 test (|Krzanowskilll990r ) to test the null hypoth- 
esis that the actual values of F2 as a function of FF, and 
the predictions of the simulated data belong to the same 
statistical sample. The T 2 statistic is defined as 



(Xi - X 2 ) 



(4) 



where Xi and X2 are the vectors containing the sample 
means and N\ — 7 and N2 = 18 are the number of real- 
izations for the data and the simulation respectively. S p is 
the pooled covariance matrix defined as S v = m§i±H2§2 
where m = Ni — 1, 712 = N2 — 1 and Si, and S 2 are the 
variance-covariance matrix for the data and model respec- 
tively. The null hypothesis is rejected for an observed value 
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Filling Factor (FF) 

Figure 2. This shows the Average Filamentarity (-F2) as a func- 
tion of Filling Factor (FF) for every alternate luminosity bins 
(Table 1). The actual data along with the corresponding results 
from ACDM N-body simulations are shown for a few representa- 
tive values of the bias. We use only the range 0.05 < FF < 0.65 
in our analysis. 
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Figure 3. This shows the measured relative bias -p- of the SDSS 
galaxies as a function of -jy . The solid curve shows the best fit 
curve to our data -pr = 0.833 + 0. 171 -jy. Fits to data from earlier 
works are also shown for comparison. The error-bars shown in the 
figure correspond to 95% confidence intervals. 

To > ( „^ + 2 "- P +i) F(p,nx +n 2 -p + 1, 1 - a) where p = 7 
is the number of variables and F is the F distribution. We 
choose a = 0.05 which allows us to reject bias values at the 
95% confidence level. For each luminosity bin, bias values 
for which the observed To is greater than T 2 — 24.758 are 
rejected with 95% confidence. 

We have tested for possible instabilities in inverting 
the pooled covariance matrix S p arising from the lim- 
ited number of realisations used to estimate Si and S2 
l|Hartlap. Simon fc Schneider! [200r3 ) by repeating the anal- 
ysis using a smaller number of realisations (Ni = 4 and 
N2 — 10). The results are found to be unchanged, except 
for a somewhat larger b range that is accepted with 95% 
confidence. 



3 RESULTS AND CONCLUSIONS 

We show results for the Average Filamentarity F2 as a func- 
tion of the Filling Factor FF in Figure [2] Increasing the 
coarse-graining causes adjacent clusters to connect up into 
progressively longer interconnected filaments and as a conse- 



Table 2. For each luminosity bin this shows the median magni- 
tude and the corresponding bias determined from our analysis. 
The errors in the bias indicate the 95% confidence interval. 



Bin 


Median magnitude 


Bias 


Bin 1 


-18.96 


0.95 ± 0.1 


Bin 2 


-19.45 


1.0 ±0.15 


Bin 3 


-20.38 


1.1 ±0.05 


Bin 4 


-21.26 


1.3 ±0.15 


Bin 5 


-21.7 


1.52 ±0.35 



quence F2 increases with FF in the simulations as well as the 
actual data. At FF in the range 0.5 — 0.6 there is a percola- 
tion transition where a large fraction of the clusters get inter- 
connected into a sing le cluster which spans across nearly the 
entire survey region dPandev fc Bharadwail[2005h . Continu- 
ing the coarse-graining considerably beyond this causes the 
clusters to thicken and the filamentary patterns are washed 
away. We restrict our analysis to the range FF < 0.65 which 
is slightly beyond the percolation threshold. A point to note 
here is that the F2 curves are sensitiv e to the area, thickness 
and number density of the sample (|Pandev fc Bharadwail 
l200rj ) which varies considerably across the luminosity bins 
(Table 1), and in the present analysis it is not meaningful to 
compare the F2 values across the luminosity bins. Consider- 
ing any one of the luminosity bins in Figure [3] we find that 
increasing the bias causes F2 to go up at low FF whereas 
the trend is reversed at large FF. This can possibly be at- 
tributed to the fact that with increasing bias the galaxies 
get concentrated into more compact regions at the expense 
of some of the large-scale coherent structures. We also note 
that the bias dependence of the F2 cur ves is somewhat differ- 
ent fr om the luminosity dependence (|Pandev fc B haradwaj 
2006) which shows a monotonic decrease with increasing lu- 
minosity. For each luminosity bin we find that the simulated 
data-set is consistent with the actual data for only a small 
range of the bias parameter. 

The results of the quantitative comparison of each lu- 
minosity bin against the biased simulations are shown in Ta- 
ble 2. It is customary to use L/L* = iq-oa(m-m') to quan- 
tify the luminosity of each bin, where M is the median mag- 
nitude of the bin an d M* has a value M* = -20.44 ± 0.01 
l|Blanton et alj 120031 ) for the SDSS. The M* galaxies lie in 
Bin 3 which has a bias value b* = 1.1. The luminosity depen- 
dence (L/L* dependence) of the relative bias b/b* is shown 
in Figure El Fitting a straight line b/b* — A + B (L/L*) by 
minimizing 

2 / 1 [(V^Data - (fo/6*)Modcl] 2 

X/V =N^ ^ (5) 

we obtain the best fit values A = 0.833 ± 0.009 and 
B = 0.171 ± 0.009 which gives an acceptable fit with 
X 1 l v = 0.066. Here we have assumed that at each luminosity 
bin fJData is half the corresponding 95% confidence interval. 
Note that such a low value of \ 2 l v i s unlikely (probabil- 
ity < 2.5%) if the errors in the different luminosity bin are 
uncorrelated and have a normal distribution. It is possible 
that these assumptions do not hold, and a further possibility 
is that using half the 95% confidence interval overestimates 

CData- 
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We find that our best fit l uminosity bias relati on is 
very close to those p roposed by iNorberg et al.l (120011 ) and 
iTegmark etHI (2004), both of which provide acceptable fits 
to our data with \ 2 l v = 0.189 and 1 .154 respective l y. The 
luminosity bias rela tion obtained by Peacock et alj (120011 ) 
through fitting the iBenoist et all (| 19961 ) data is consider- 
ably steeper. Their luminosity bias relation shows consid- 
erable deviations from our data at both the low and high 
luminosity end s and is ruled out wit h y 2 jv = 4.697. This is 
consistent with lNorberg et alj (|200ll ) who find considerable 
deviations between their data and the IBenoist et al l (|l996l ) 
data at high luminosities. 

The observed luminosity bias relation is an useful in- 
put for models of galaxy formation. Hiera rchical models 
of galaxy formation (e.g. IWhite et al.l Il987h generally as- 
sume the brighter galaxies to be associated with the more 
massive dark matter halos which have a stronger clus- 
tering. The s emi-a nalytic model of galaxy formation of 
iBenson et al. is consistent with the observed luminos- 

ity bias relation IjNorberg et al. I l200ll ). IZehavi et al.l (|2005l ) 
show that the halo model with suitably chosen parameters 
can correctly predict the observed luminosity dependence. 
We note that the focus of much of these works has been 
on the luminosity dependence of the two-point statistics. 
The largest currently available implementation of a semi- 
analytic model of g alaxy formation (The Millennium Run, 
ISpringel et al.l l2005) fails to correctly predict the luminosity 
dependence of the filamenta rity for the brightest galaxies 
l|Pandev fc Bharadwaill2006l ). Finally, it is interesting that 
our analysis which uses a global property of the galaxy distri- 
bution, namely its morphology, predicts the same luminosity 
bias relation as determined from the two point statistics. 
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